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The Study of Stationarity of One-dimensional Stochastic Processes 


Modele ARIMA (p,d,q) są obecnie szeroko stosowane do opisu dyna- 
miki prognozowania procesów losowych, pomimo że metody identyfi- 
kacji! tych modeli są niejednoznaczne i często wysoce wątpliwe. Podsta- 
wę identyfikacji wspomnianych modeli stanowią uwagi zawarte w pracy 
G. E. P. Boxa i G. M. Jenkinsa [6, s. 174—187], zgodnie z którymi ocena 
parametrów p, d, q może być dokonana w oparciu o analizę wykresów 
ocen wartości funkcji autokorelacji i autokorelacji cząstkowej. W niniej- 
szej pracy uwagę skoncentrowano na problemie doboru rzędu różnicowa- 
nia d, czyli na badaniu stacjonarności w szerszym sensie procesu repre- 
zentowanego przez określony szereg czasowy. 

Praca składa się z sześciu części i zakończenia. Po uwagach wstępnych, 
w drugiej części pracy przedstawiono uwagi dotyczące estymacji pod- 
stawowych charakterystyk procesów losowych. W części trzeciej omówio- 
no istotę zaproponowanej przez Boxa i Jenkinsa metody badania stacjonar- 
ności. Część czwartą poświęcono wykorzystaniu medianowego testu serii 
do badania stacjonarności jednowymiarowych ciągów losowych. W części 
piątej dokonano oceny efektywności metod badania stacjonarności, które 
są stosunkowo często wykorzystywane w praktyce, czyli metody Boxa 
i Jenkinsa oraz medianowego testu serii. Ocenę tę przeprowadzono w opar- 
ciu o wyniki badań symulacyjnych. W szóstej części pracy zapropono- 
wano metodę badania stacjonarności, która — jak sądzi autor — umożli- 
wia bardziej wiarygodną ocenę stacjonarności, niż metody dotychczas 
stosowane. 


1 Identyfikacja modeli ARIMA (p, d, a) polega na ocenie parametrów p, d i a. 
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ESTYMACJA CHARAKTERYSTYK SŁABO STACJONARNYCH 
PROCESÓW LOSOWYCH 


Pod pojęciem słabo stacjonarnego procesu losowego rozumiemy pro- 
ces stacjonarny ze względu na wartość oczekiwaną i funkcję kowariancji, 
czyli: 

E(X(t)) =p (1) 
A 
t=0,+1,+2,... 
A |t,— te] =|ty—t,] (= eov( X(t) X(te) = cov(X (ts), X(t,)). 
tstatatı=0, +1, +2,... (2) 

W przypadku stacjonarnosci procesu ze wzgledu na funkcje kowa- 
riancji mamy możliwość * wprowadzenia pojęcia funkcji autokowariancji 


o postaci: 
V(t) =cov(X(t.),K(ty)), (3) 
gdzie: t=t— ty(r=0, + 1, + 2,...) 


Z własności (2) wynika parzystość funkcji autokowariancji, czyli 
v(r)"Y(— 1). 

Przebieg procesów losowych charakteryzowany jest często na podsta- 
wie ich funkcji korelacji, czyli 

e(t,t;) =cov(X(t,),X(t;))/(s(t,) -s(t,)), (4) 
gdzie: 
s(tt)= 4/ cov(X(t,),X(t)) 

Proces losowy jest stacjonarny ze wzgledu na funkcje korelacji, gdy 
funkcję o(ti,t;) można zastąpić funkcją jednej zmiennej 1=t,—ti(t,,tj=0, 
+1, +2,..). Funkcję o(t) nazywamy funkcją autorelacji. Przedstawione 
uwagi wskazują, że stacjonarność ze względu na funkcję korelacji nie im- 
plikuje stacjonarności ze względu na kowariancje *. 

Omówione wyżej charakterystyki szacowane są w przypadku założe- 
nia o słabej stacjonarności i ergodyczności danego procesu losowego na 
podstawie następujących wzorów: 


x=n' VXT (5) 


t=l 


2 W pracy świadomie użyto słów „mamy możliwość”. W wielu pracach z zakre- 
su analizy szeregów czasowych pojęcie funkcji autokorelacji jest bowiem używane 
dla procesów niestacjonarnych ze względu na funkcję kowariancji. Tego typu po- 
dejście jest — według autora niniejszej pracy — niezgodne z teorią procesów lo- 
sowych. 

3 Uwaga ta wynika z faktu, że można wyobrazić sobie taki proces losowy, dla 
którego o (ti, ty) nie zależy od tı, ty, lecz od t=tj—tj, pomimo, iż proces ten nie 
jest stacjonarny ze względu na funkcję kowariancji. 


Badanie stacjonarności jednowymiarowych procesów... 423 


RE: 


e= ) Dn), (6) 
t=1 
T(t) =c(z)/c(0), t=0,1,....m 
1=0,1,..,m (7) 
gdzie: (x,Xa,..,Xy) Oznacza pojedynczą realizację analizowanego procesu 


losowego, 
m jest punktem odcięcia funkcji autokowariancji i autokorelacji. 
Średni błąd szacunku parametru g(r) oceniany jest w oparciu o for- 
mułę Bartletta [6, s. 179]. Formuła ta przedstawiona jest w następujący 
sposób: 
s(r(k)) =n"*(1+ 2(r*(1)+ r*(2)+...+ r*(q)))*2, 
k>q (8) 
Wartości funkcji autokoleracji uznawane są jako statystycznie nie- 
istotne, gdy 
r(k)EL—2-s(r(k)),2*s(r(k))). (9) 


OCENA RZĘDU RÓŻNICOWANIA W METODZIE BOXA-JENKINSA 


Załóżmy, ze do szeregu czasowego X,,X2,....X, stanowiącego realizację 
procesu losowego (X(t); 
t=0, +1,+2,...) dopasowujemy model ARIMA, p, d, q) o postaci: 
F(B)V °x, "Q(B)a+, (10) 
gdzie: 
F(B)=l—ayB —ap: B*—...— Opp BP, 
Q(B)F1—haB —ha B?—...—haa BY, 
V°x,=(1—B)? x, 
Bxt= Xt, Bar, 
BPxt= Xp  Biap=ary, 
{at} oznacza proces czysto losowy o wariancji ôa. 

Identyfikacja modeli (10) polega na ocenie parametrów p, d, q ozna- 
czających odpowiednio: rząd składowej autoregresyjnej (AR), rząd róż- 
nicowania doprowadzającego proces (X+) do słabej stacjonarności i rząd 
składowej występującej w postaci średniej ruchomej (MA). 

G. E. P. Box i G. M. Jenkins proponują, by rząd różnicowania oznaczał 
taką najmniejszą liczbę z ciągu 0, 1, 2,.., dla której szereg czasowy 
w+= Vdxętmożna uznać za realizację słabo stacjonarnego procesu losowe- 
go. Wspomniani autorzy sugerują, że decyzja odnośnie stacjonarności 
(bądź niestacjonarności) danego procesu losowego może być podejmowana 
w oparciu o obserwację ,,zachowania” ocen wartości funkcji autokore- 
lacji danego procesu. Powyższa sugestia wynika z następującego stwier- 
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dzenia [6, s. 175]: „w przypadku modelu stacjonarnego, w którym żaden 
z pierwiastków nie leży blisko okręgu jednostkowego funkcja autokore- 
lacji szybko zanika dla dużych i średnich k”. Stąd wyciągany jest wnio- 
sek, że w sytuacji, gdy funkcja autokorelacji nie wykazuje tendencji do 
szybkiego zaniku analizowany szereg czasowy można rozważać jako reali- 
zację procesu niestacjonarnego. Box i Jenkins podejmują próbę formal- 
nego wyjaśnienia przedstawionych wyżej spostrzeżeń. Podstawę tego wy- 
jaśnienia stanowi fakt, że funkcja autokorelacji słabo stacjonarnego pro- 
cesu losowego ARMA (p, q) spełnia równanie różnicowe o postaci: 

o(B)o(k)=0 k>q, (11) 
gdzie 


9(B) = II(1—G,B). 
i =1 


Przy założeniu, że wszystkie pierwiastki równania różnicowego (11) 
są jednokrotne, jego rozwiązanie może być przedstawione w następujący 
sposób: 

o(k) =A,G¥,+ AG", +... + A,G*, k>q—p, (12) 

Z równania (12) wynika, że w przypadku, gdy przynajmniej jeden 
pierwiastek równania charakterystycznego leży w pobliżu okręgu jedno- 
stkowego, funkcja autokorelacji opada wolno ([6, s. 176]). Uważam, ze 
równanie (12) nie może być wykorzystywane do wyciągania jakichkol- 
wiek wniosków odnośnie przebiegu teoretycznej funkcji autokorelacji 
procesu niestacjonarnego, gdyż funkcja taka nie istnieje. Można jedynie 
podjąć próbę analizy kształtowania się wartości estymatora funkcji auto- 
korelacji procesu stacjonarnego w sytuacji, gdy dysponujemy realizacją 
procesu niestacjonarnego. Wszelkie uogólnienia mogą być dokonywane 
jedynie na podstawie wyników badań symulacyjnych. 


MEDIANOWY TEST SERII W BADANIU STACJONARNOŚCI 
JEDNOWYMIAROWYCH CIĄGÓW LOSOWYCH. 


Obok metody Boxa i Jenkinsa podstawę badania stacjonarności pro- 
cesów losowych stanowi również medianowy test serii omówiony m. in. 
w [4] i [7]. Zgodnie z propozycją J.S. Bendata i A. G.Piersola badanie 
stacjonarności sprowadzane jest do oceny stacjonarności ze względu na 
wartość średnią i wariancję. Zastosowanie wspomnianego wyżej testu do 
weryfikacji hipotezy o stacjonarności ze względu na średnią wymaga po- 
działu n elementowego szeregu czasowego (realizacja procesu) na k roz- 
łącznych s elementowych podciągów (x, Xz. Xs}, {Xs¢1 Xstayey Kas), 
{Xa n-s+1r--+Kk.s |. Dla każdego podciągu wyznaczane są wartości 
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, 
zj =s- DZE (i=1, 2,..,k), (13) 
j=l 

Przy założeniu prawdziwości H, o stacjonarności ze względu na śred- 
nią, wartości Xx; i=1, 2,..,k stanowią realizacje zmiennej losowej o war- 
tości oczekiwanej p. Widać więc, że ciąg wartości X; i=1, 2,..,k powinien 
charakteryzować się losowością. Stąd liczba serii w ciągu uzyskanym 
w wyniku porównania wartości x, z ich medianą (Me,) powinna być zbli- 
żona do liczby serii charakterystycznej dla ciągu niezależnych obserwacji 
zmiennej losowej. Sprawdzianem w omawianym teście jest liczba serii 
uzyskanych w wyniku porównywania x, z Me,. Przy budowie obszaru 
krytycznego wykorzystywane są następujące równości: 

P {1<l}= 4a i P{1<1,}=1— ha, 
gdzie: 
l oznacza liczbę serii uzyskanych w danym ciągu, l, l; oznaczają 
wartości krytyczne odczytane z tablic rozkładu liczby serii dla po- 
ziomu istotności a (zob. [7]). 
W świetle powyższych uwag widać, że gdy wyznaczona wartość 1 spełnia 
jedną z nierówności 1<], lub 1>1, to hipotezę o stacjonarności rozważa- 
nego procesu należy odrzucić. W przypadku gdy 1,<|<1,, nie mamy pod- 
staw do odrzucenia H,. 

Badanie stacjonarności procesu ze względu na wariancję wymaga ana- 
logicznego do przedstawionego wyżej postępowania z ciągiem wariancji 
s4; i=1, 2,..,k. Proces uznajemy za stacjonarny, gdy nie mamy podstaw 
do odrzucenia hipotezy H, o jego stacjonarności zarówno ze względu na 
średnią, jak i na wariancję. 


OCENA EFEKTYWNOŚCI KLASYCZNYCH METOD BADANIA 
STACJONARNOSCI JEDNOWYMIAROWYCH CIAGOW LOSOWYCH 


W tej części pracy oceniona zostanie AUG ank mediano- 
wego testu serii i metody Boxa-Jenkinsa do badania stacjonarności proce- 
sów losowych. Statystyczna ocena wymienionych metod zostanie doko- 
nana w oparciu o wyniki badań symulacyjnych. 

Ocena obydwu sposobów postępowania wymagała generowania reali- 
zacji wybranych słabo stacjonatnych i niestacjonarnych procesów AR(p). 
Rese te ORA w EL zastosowania następującej formuły: 

in 
t=ptl,....n 
2k t=1, 2,.., p (14) 


SERES 
Tew, 
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gdzie: 
api; 1=1, 2,...,.p — parametry autoregresji, 
(e; t=pt1,...n) — liczby pseudolosowe o rozkładzie normalnym 
N(0, 1). 


Liczby pseudolosowe o rozkładzie normalnym N (0,1) uzyskiwano 
w wyniku zastosowania generatora wykorzystującego centralne twier- 
dzenie graniczne. Opis tego generatora można znaleźć np. w [9, s. 84]. 
W przypadku, gdy ocenie poddawano metodę Boxa-Jenkinsa generowano 
realizacje następujących procesów: 


AR(1):a1u=0,3;  au=0,4; au=0,9; au=1l,0;  a1=1,05 


AR(2)‘an=1, a,,.=—0,9 — proces stacjonarny 
aa =1,1;  a.=—0,1 — proces niestacjonarny 
AR(3)!an=0,8;  a3+=0,51;  a;,=— 0,378 — proces stacjonarny 
@3,=0,9; «3470,52;  a3—0,42 — proces niestacjonarny 


Pojedynczy eksperyment polegał na 100-krotnym generowaniu reali- 
zacji jednego z wyżej wymienionych procesów. Na podstawie każdej rea- 
lizacji dokonywano oceny wartości funkcji autokorelacji rzędu 1, 2,...,20 
i okleślano rząd począwszy, od którego wartości funkcji autokorelacji 
okazały się statystycznie nieistotne. Rząd ten określany był na podstawie 
wyrażenia (9). W oparciu o każde 100 wyników wyznaczano średni rząd 
począwszy od którego wartości funkcji autokorelacji okazywały się sta- 
tystycznie nieistotne. 


Średni rząd obliczono na podstawie następującego wzoru: 


i=1 
sr=1007 Sik, (15) 
100 


gdzie: l, oznacza rząd począwszy od którego wartości funkcji autokorelacji 
okazały się nieistotne statystycznie w i-tym powtórzeniu. 


Odchylenie standardowe wyników uzyskanych w danym eksperymen- 
cie obliczano według wzoru o postaci 


100 
1007: X (Li —sr)ż. (16) 
Lm 


i=l 


Dla kadżego procesu generowano realizacje o liczbie obserwacji n=50, 
100, 200 i 500. 


Wyniki badań przeprowadzonych w omówiony wyżej sposób zawarto 
w tabelach 1 i 2. 


Analiza uzyskanych wyników pozwala zauważyć, że o ile w przypad- 
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Tab. 1. Ocena metody Boxa-Jenkinsa 
The estimation of Box-Jenkins’ method 


Liczba obserwacji 


50 100 200 500 
Nazwa 7 ——— I, aa |b ———_ ZAK Fs 
procesu odchy- odchy- odchy- odchy- 
średni lenie średni lenie średni lenie średni lenie 
rząd |standar-| rząd  |standar-| rząd standar-| rząd standar- 
dowe dowe | | dowe | dowe 
AR (1) 
an=0,3 3,63 2,76 5,53 5,61 6,08 5,68 8,36 6,59 
AR (1) 
a11=0,4 3,42 2,92 4,21 4,01 6,26 5,64 6,36 5,53 
AR (1) 
an =0,5 3,39 2,82 4,42 4,02 5,37 4,60 5,46 4,99 
AR (1) 
a11=0,6 3,23 4,15 5,19 4,15 5,29 4,08 7,12 5,00 
AR (1) 
a11=0,7 3,59 2,94 5,19 3,66 5,73 3,46 7,46 3,53 
AR (1) 
a11=0,8 3,29 1,29 4,88 2,12 6,33 3,21 8,38 3,30 
AR (1) 
a1=0,9 3,69 1,62 6,13 2,42 8,96 2,52 12,52 3,01 
AR (1) 
au=1 4,31 0,64 8,03 1,33 14,61 3,42 — — 
AR (1) 
«11=1,05 4,79 0,52 8,01 0,36 — — — = 


Żródło: Obliczenia własne. 


c — wszystkie oceny okazały się statystycznie istotne. 


ku prób bardzo dużych (500 obserwacji) istnieje możliwość jednoznacz- 
nego stwierdzenia czy dany proces można uznać za stacjonarny, czy nie‘, 
to w przypadku prób o liczbie obserwacji n=50 i n=100 wyciągnięcie ta- 
kiego wniosku jest praktycznie niemożliwe. Powyższa uwaga wynika 
z tego, że dla prób o liczbie obserwacji 50 i 100, średni rząd powyżej któ- 
rego wartości funkcji autokorelacji są statystycznie nieistotne jest zbliżo- 
ny zarówno dla procesów stacjonarnych, jak i niestacjonarnych. Jest on 


4 W przypadku, gdy generowano bardzo długie szeregi czasowe stanowiące 
realizacje procesów niestacjonarnych, to wszystkie wartości funkcji autokorelacji 
okazywały się statystycznie istotne. 


28* 
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wprawdzie dla procesów niestacjonarnych wyższy lecz tak nieznacznie, 
iż trudno jest wyciągnąć tutaj wnioski, które mogłyby być przydatne 
w praktyce badania stacjonarności. Uzyskane wyniki uświadamiają w peł- 
ni, jak bardzo zawodne jest wykorzystywanie metody Boxa i Jenkinsa 
do badania stacjonarności. 

Kolejno przystąpiono do oceny efektywności badania stacjonarności 
przy wykorzystaniu medianowego testu serii. Celem badań było określe- 
nie mocy medianowego testu serii* w sytuacji, gdy stosowany jest on do 
badania stacjonarności procesu losowego. W tym przypadku zgodnie 
z zasadami podanymi na wstępie niniejszej części pracy generowano rea- 
lizacje niestacjonarnych procesów losowych AR(1) z ay =1 i ay =1,1. Po- 
jedynczy eksperyment polegał na 100 krotnym wygenerowaniu realizacji 


Tab. 2. Ocena metody Boxa-Jenkinsa 
The estimation of Box-Jenkins’ method 


Liczba obserwacji 


50 100 200 500 
Nazwa | === o | 6 eel pcs, 
procesu | odchy- odchy- odchy- odchy- 
średni | lenie | średni | lenie | średni | lenie | średni | lenie 
„ rząd |standar-| rząd .|standar-| rząd |standar-| rząd jstandar- 
dowe dowe dowe dowe 
AR (2) 
a1=1 3,86 1,78 6,33 2,75 8,48 2,34 12,68 2,74 
az= —0,09 ; 
AR (2) 
an=1,1 4,29 0,97 7,69 159 : 145 2,90 — = 
a22=—0,1 
AR (3) 
an=0,8 
as2=0,51 3,93 1,00 6,31 1,50 9,39 2,74 118 3,72 
as = —0,378 e è: Ah A 
AR (3) l l 
an =0,9 
as=0,52 4,43 0,75 7,99 1,30 ..: 15,18 3,21 He — 
asa — 0, 42. i 


c — wszystkie oceny reas sie statystycznie istotne. 
Źródło: obliczenia własne. 


5 Pod pojęciem mocy. testu rozumie się tu, prawdopodobieństwo, odrzucenia hi- 
potezy głoszącej stacjonarność procesu losowego w. przypadku, gdy realizacja re- 
prezentowała proces niestacjonarny. 
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Tab. 3. Moc testu stacjonarności 
The power of the test of stationarity 


Nazwa Moc Nazwa Moc 
procesu testu procesu testu 
AR (1) N=50 AR (1) N=50 
an=1 M=5 0,874 a=11 M=5 0,856 

N=50 
M=10 0,79 
N=100 N=100 
M=5 0,972 M=5 0,986 
N=100 N=100 
M=3 0,858 M=8 0,846 
N=100 N=100 
M=10 0,850 M=10 0,844 
N=150 N=150 
M=5 1 M=10 0,99 
N=150 N=150 
M=10 l M=15 0,854 
N=150 
M=15 0,892 
N=200 
M=20 0,872 


Żródło: obliczenia własne. 


wymienionych wyżej procesów. Badanie stacjonarności przeprowadzano 
dla szeregów czasowych o różnej długości (n=50, 100, 150, 200) i przy 
założeniu różnej długości odcinków, na które szeregi te były dzielone. 
Wyniki badań przedstawiono w tabeli 3. 

Rezultaty, które uzyskano w wyniku badania mocy testu stacjonar- 
ności w przypadku realizacji niestacjonarnych procesów AR(1) wydawały 
się być obiecujące. We wszystkich bowiem przypadkach prawdopodobień- 
stwo odrzucenia fałszywej hipotezy zerowej było wysokie (0,85 — 1). Po- 
nieważ przy postępowaniu polegającym na dwukrotnym zastosowaniu 
testu serii nie ma możliwości określenia prawdopodobieństwa popełnienia 
błędu I rodzaju, zdecydowano się na podjęcie badań symulacyjnych, które 
umożliwiłyby oszacowanie tego prawdopodobieństwa. Badania ograniczo- 
no jedynie do realizacji procesu AR(1) z a, =0,8, gdyż już pierwsze wy- 
niki przekreśliły praktyczną przydatność analizowanego testu. Okazało 
się bowiem, ze dla n=50 i m=5 stacjonarność została stwierdzona jedy- 
nie w 21 na 100 powtórzeń, dla n=100, m=10 w 22 powtórzeniach, dla 
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Tab. 4. Empiryczny rozkład estymatora parametru au 
Empirical distribution of the estimator of parameter an 


AR(1) a11=—0,99 AR (1) au=—1 AR (1) an= —1,01 


Przedział — 

klasowy n=50 n=100 n=200 |n=50 n=100 n=200 |n=50 n=100 n=200 

| n=300 n=300 n=300 

poniżej (—1,175) o 0 0 0 0 0 0 0 0 0 0 9 
(—1,175)—( — 1,15) 0 0 0 © 0 0 0 0 0 0 0 0 
(—1,15——1,125) 0 0 0 0 0 0 0 0 0 0 0 0 
(—1,125—(—1,1) 0 0 0 0 0 0 0 0 0 0 0 D 
(—1,1)4—1,075) 0 0 0 0 0 0 0 0 0 0 0 0 
(—1,075)—{ — 1,05) 0 0 0 0 0 0 0 0 0 0 0 0 
(—1,05){— 1,025) 3 0 0 0 2 0 0 0 3 0 0 0 
(—1,025—(— 1,00) 4 10 4 1 35 34 37 27 74 80 93 93 
(—1,00——0,975) 30 56 71 82 38 53 56 73 13 11 6 6 
(—0,975)—(—0,95) 24 2% 18 16 11 6 7 0 6 5 1 1 
(—0,955— — 0,925) 9 13 6 1 8 2 0 0 2 1 0 0 
(—0,925)— — 0,9) 4 1 1 0 2 4 0 0 0 3 0 0 
(—0,9)—{— 0,875) 1 0 0 0 1 1 0 0 1 0 0 0 
(—0,875)—{ — 0,85) 1 0 0 0 1 0 0 0 0 0 0 0 
(—0,85)—{ — 0,825) 0 0 0 0 1 0 0 0 0 0 0 0 
(—0,825)—{ — 0,8) 2 0 0 0 1 0 0 0 1 0 0 0 
powyzej (—0,8) 2 0 0 0 0 0 0 0 0 0 0 0 


a liczba ta określa liczbę powtórzeń, w których estymator parametru au przyjął 
wartość poniżej (—1,175). 


Źródło: Obliczenia własne. 


n=l100, m=15 w 12 powtórzeniach, zaś dla n=300, m=15 zaledwie w 1 
powtórzeniu. W świetle powyższych uwag możemy stwierdzić, że stoso- 
wane dotychczas metody badania stacjonarności procesów losowych są 
wysoce zawodne. O ile jednak metoda Boxa-Jenkinsa ma walory prak- 
tyczne w przypadku szeregów bardzo długich, to omówiony w pracy test 
stacjonarności jest pozbawiony takich walorów. 


PROPOZYCJA METODY BADANIA STACJONARNOŚCI 
PROCESÓW LOSOWYCH 


Podstawę proponowanej w nijnieszej pracy metody badania stacjonar- 
ności stanowi hipoteza, że w przypadku rozważania niestacjonarnego 
procesu losowego jako procesu AR(1), uzyskana klasyczną metodą naj- 
mniejszych kwadratów ocena parametru a, przyjmuje wartość większą 
lub równą 1. Powyższą hipotezę weryfikowano w oparciu o wyniki badań 
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Tab. 5. Empiryczny rozkład estymatora parametru au 
Empirical distribution of the estimator of parameter au 


AR (1) a11=0,99 AR (1) au=1 | AR (1) a1=1,01 

Przedział a | mm 2. 
klasowy n=50 n=100 n=200 |n=50 n=100 n=200 |n=50 n=100 n=200 
n=300 n=300 n=300 
poniżej 0,8 12 0 0 0 7 1 0 0 9 0 0 0 
0,8—0,825 7 0 1 0 4 1 0.0 4 0 0 0 
0,825—0,85 6 2 0 0 10 4 0 o 2 2 0 0 
0,85—0,875 7 5 0 0 12 4 1 0 3 1 0 0 
0,875—0,9 11 3 3 1 8 4 1 0 8 0 1 0 
0,9—0,925 6 9 2 0 1 2 1 0 15 1 0 0 
0,925—0,95 22 27 12 14 #2 20 6 6 9 8 0 0 
0,95—0,975 16 28 37 34 14 23 31 21 15 15 5 3 
0,975—1 10 24 44 49 8 27 55 64 17 31 10 7 
1—1,025 3 2 1 3 3 4 5 9 15 42 84 = 90 
1,025—1,05 0.0 0 0 0.0 0 0 2 0 0 0 
1,05—1,075 0 0 0 0 0.0 0 0 1 0 0 0 
1,075—1,1 0 0.0 0 0 0.0 0 0.0 0 0 
1,1—1,125 0.0 0 0 0.0 0 0 0 0 0 0 
1,125—1,15 0.0 0 0 0 0.0 0 0 0 0 0 
1,15—1,175 0 0 0 0 0.0 0 0 0.060 0 0 


Źródło: Obliczenia własne. 


symulacyjnych przeprowadzonych dla wybranych stacjonarnych i niesta- 
cjonarnych autoregresyjnych procesów losowych. 

Procesy te dobierano w ten sposób, by pierwiastki odpowiadających 
im równań charakterystycznych znajdowały się w pobliżu lub na granicy 
stacjonarności. Zgodnie z formułą (14) generowano realizacje następują- 
cych procesów losowych: 


AR(1):au= 0,99; n=l; a,,=1,01 
au=— 0,99; a,=—1; a,=—1,01 
AR(2):a4=1,1; a4=—0,1 — proces niestacjonarny (pierwiastek na 


granicy stacjonarnosci), 
@,=1,1; a= —0,09 — proces niestacjonarny, 
Q,=1,1; @=—0,11 — proces stacjonarny, 
aan =— 0,21;  ay+=0,808 — proces niestacjonarny, 
ax: — 0,2; dz =0,8  — proces niestacjonarny (pierwiastek 
na granicy stacjonarności 
an =— 0,19;  a„=0,792 — proces stacjonarny 
AR(3):a4=1,1; 4450,32;  a;,=—0,42 — proces niestacjonarny 
(pierwiastek na granicy 
stacjonarności) 
ay=1,11; 4,.=0,319; a„=— 0,4242 — proces niestacjonarny 
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Tab, 6. Empiryczny rozkład estymatora parametru au 
Empirical distribution of the estimator of parameter au 


AR (2) asi=1,1 AR (2) asi==1,1 AR (1) au=1,1 
Przedział az= — 0,11 az= —0,1 a22= —0,09 

klasowy n=50 n=100 n=200 |n=50 n=100 n=200 |n=50 n=100 n=200 

n=300 n=300 n=300 
poniżej 0,8 1 0 0 0 4 0 0 0 8 0 0 0 
0,8—0,825 2 1 0 0 4 1 0 0 1 0 0 0 
0,825—0,85 8 0 0 0 4 1 0 0 0 0 0 0 
0,85—0,875 4 4 0 0 2 3 0 0 8 0 0 0 
0,875—0,9 14 4 0 0 138 3 0 0 6 2 0 0 
0,9—0,925 14 13 3 1 0 9 2 0 8 4 0 0 
0,925—0,95 18 25 5 3 15 11 6 0 11 1 1 0 
0,95—0,975 15 29 30 28 16 33 27 2% #17 9 4 0 
0,975—1,0 9 22 61 68 15 37 #57 79 23 22 7 4 
1,0—1,025 5 2 1 0 7 2 8 1 20 61 88 96 
1,025—1,05 0.0 0 0 0 0 6 0 3 1 0 0 
1,05—1,075 o o o 0 0 0 0 0 0 6 0 © 
1,075—1,1 0.0.0 0 0 0 0 0 0 0 0 0 
1,1—1,125 o o o o 0 0 0 0 0 0 0 © 
1,125—1,15 0.0.0 0 0 0 o 6 0 0 0 d 
1,15—1,175 0.0 o o 0 0 0 0 0 0 0 0 


Źródło: Obliczenia własne. 


a,=1,11;  4+=0,319; a; —0,4242 — proces niestacjonarny 
ay =l; a= 0,33; Qa; =— 0,37. — proces stacjonarny 

Pojedynczy eksperyment polegał na 100-krotnym wygenerowaniu rea- 
lizacji jednego z wymienionych wyżej procesów. Na podstawie każdej 
realizacji w oparciu o metodę najmniejszych kwadratów szacowano pa- 
rametr a, modelu AR(1). Uwzględniając wyniki dla wszystkich 100 po- 
wtórzeń uzyskiwano empiryczne rozkłady estymatora parametru oy. Ba- 
dania przeprowadzono dla prób o liczebnościach n=50, 100,,200 i 300. 
Uzyskiwane wyniki zamieszczono w tabelach 4, 5, 6, 7 18. 

Szczegółowa analiza informacji zawartych w tabelach 4—8 pozwala 
stwierdzić, że dla 200 i 300 elementowych realizacji procesów niestacjonar- 
nych o pierwiastkach równania charakterystycznego znajdujących się wew- 
nątrz okręgu o promieniu jednostkowym w. około 90°/o eksperymentów 
uzyskiwano ocenę a, powyżej wartości 1 lub ponieżj wartości — 1, Dla 
prób 50 i 100 elementowych wyniki te nie były już tak korzystne i kształto- 
wały się odpowiednio na poziomie około 15—-20°/o oraz 42—60°/o. Widać jed- 
nak wyraźnie, że w przypadku, gdy bezwzględna wartość z oceny parame- 
tru a, jest większa lub równa 1, to proces reprezentowany przez daną pró- 
bę można zakwalifikować do klasy niestacjonarnych. Rozkład estymatora 
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Tab. 7. Empiryczny rozkład estymatora parametru au 
Empirical distribution of the estimator of parameter an 


AR (2) az1= —0,21 AR (2) az1=—0,2 AR (2) as1=—0,19 


Przedziat a22= 0,808 Q22=0,8 a22=0,792 

klasowy n=50 n=100 n=200 |n=50 n=100 n=200 |n=50 n=100 n=200 
n=300 n=300 n=300 

poniżej (—1,175) 0.0 0 0 0 0 0 0 0 0 0 0 
(—1,175)——1,15) 0 0 0 0 0 0 0 0 0 0 0 0 
(--1,15— —1,125) o o 0 0 0 0 0 0 6 0 0 0 
(—1,125)—— 1,1) 0.0 0 0 0 0 0 0 0 0 0 0 
(—1,1)—(—1,075) 0.0 0 0 0 0 0 0 0 0 0 0 
(—1,075)——1,05) 0.0 0 0 0 0 0 0 0 0 0 0 
(—1,055— — 1,025) 0.0 0 0 0 0 0 0 0 0 0 0 
(—1,025)—(—1,0) 35 61 7 9 4 #7 2 3 1 0 0 0 
(—1,00—(—0,975) 386 16 10 5 23 36 40 48 6 10 2 1 
(—0,975——0,95) 7 7 2 2 19 13 21 2 7 8 10 15 
(—0,95—(—0,925) 7 3 4 2 9 7 1 7 15 15 19 16 
(—0,925)——0,9) 2 1 3 0 4 8 7 5 10 9 13 19 
(—0,9)——0,875) 2 3 3 0 2 7 2 5 9 10 10 19 
(—0,875)—(—0,85) 1 1 0 0 6 5 1 4 10 8 8 8 
(—0,85)—(— 0,825) 2 2 1 1 3 2 4 3 3 7 9 7 
(—0,825)— —0,8) 1 3 0 0 4 2 2 1 2 4 4 4 
powyżej (—0,8) 7 3 2 0 2 13 0 4 37 29 25 11 


Źródło: Obliczenia własne. 


uzyskany dla procesów stacjonarnych wyraźnie wskazuje, że przy wyko- 
rzystaniu powyższej uwagi do badania stacjonarności prawdopodobieństwo 
popełnienia błędu polegającego na zakwalifikowaniu procesu stacjonarne- 
go do klasy niestacjonarnych jest bardzo małe i należy do przedziału od 
0 do 0,05. Łatwo zauważyć, że w przypadku procesów niestacjonarnych 
dla prób o liczbie obserwacji 200 i 300 prawdopodobieństwo tego, że bez- 
względna wartość oceny parametru a, jest mniejsza od 0,95 jest małe 
i należy do przedziału od 0 do 0,1. Dla 100 prób elementowych małe jest 
prawdopodobieństwo tego, że bezwzględna wartość oceny parametru a;, jest 
mniejsza od 0,925. Należy ono do przedziału od 0 do 0,2. W przypadku 
prób 50 elementowych trudno jest mówić o jakiejkolwiek prawidłowości. 

W świetle powyższych uwag można stwierdzić, że zarówno dla prób 
małych, jak i dużych wystąpienie większej lub równej 1, bądź mniejszej 
lub równej —1 wartości estymatora parametru a, pozwala stwierdzić 
niestacjonarność procesu reprezentowanego przez dany szereg czasowy 
z niewielkim prawdopodobieństwem popełnienia błędu. W przypadku 
prób o liczbie obserwacji większej lub równej 100 wystąpienie bezwzglę- 
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Tab. 8. Empiryczny rozkład estymatora parametru an 
Empirical distribution of the estimator of parameter au 


AR (3) asi=1 AR (3) as1=1,1 AR (3) as1=1,11 
as2=0,33 ase =0,32 a32=0,32 
Wartość aas= — 0,378 ass = — 0,42 as3= —0,42 
estymatora I | E IE 
n=50 n=100 n=200 |n=50 n=100 n=200 |n=50 n=100 n=200 

n=300 n=300 n=300 
poniżej 0,8 2 0 0 0 0 0 0 0 0 0 0 0 
0,8—0,825 1 0 0 0 2 0 0 0 0 0 0 0 
0,825—0,85 1 0 0 0 1 0 0 0 3 0 0 0 
0,85—0,875 8 3 0 0 4 0 0 0 5 0 0 0 
0,875—0,9 11 3 0 0 6 0 0 0 2 0 0 0 
0,9—0,925 16 17 2 0 8 2 0 0 6 0 0 0 
0,925—0,95 19 21 21 13 21 6 0 0 12 3 0 0 
0,95—0,975 23 41 67 79 19 24 8 0 28 6 1 0 
0,975—1 14 15 10 8 21 48 66 78 24 32 1 4 
1—1,025 3 0 0 0 13 19 26 22 11 52 92 96 
1,025—1,05 1 0 0 0 4 1 0 0 8 7 0 0 
1,05—1,075 1 0 0 0 1 0 0 0} 1 0 0 0 
1,075—1,1 0 0 0 0 0 0 0 0 0 0 0 0 
1,1—1,125 0 0 0 0 0 0 0 0 0 0 0 0 
1,125—1,15 0 0 0 0 0 0 0 0 0 0 0 0 
1,15—1,175 0 0 0 0 0 0 0 0 0 0 0 0 


Żródło: Obliczenia własne. 


dnej wartości estymatora parametru a, na poziomie niższym od 0,925 
pozwala zakwalifikować proces reprezentowany przez dany szereg czaso- 
wy do klasy stacjonarnych. Przedstawione powyżej wnioski stanowią 
przekonywającą podstawę proponowanej w pracy metody badania sta- 
cjonarności. 

Metoda ta nie rozstrzyga problemu stacjonarności bądź niestacjonar- 
ności w sytuacji, gdy bezwzględna wartość oceny parametru a, należy 
do przedziału (0,925; 1). Proponuje się, by w takiej sytuacji proces był 
niestacjonarny, czyli, aby badaniu stacjonarności zgodnie z przedstawio- 
nymi wyżej uwagami poddany został szereg pierwszych bądź kolejnych 
różnic. 


ZAKOŃCZENIE 


Przeprowadzone w pracy badania upoważniają do sformułowania na- 
stępujących wniosków: 

1. Medianowy test serii nie daje możliwości wiarygodnego rozstrzy- 
gania problemu stacjonarności bądź niestacjonarności procesów losowych; 
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2. Rozróżnienie realizacji procesów losowych stacjonarnych i niesta- 
cjonarnych w oparciu o analizę zachowania wartości estymatora funkcji 
autokorelacji jest w przypadku realizacji 50, 100 i często nawet 200 ele- 
mentowych praktycznie niemożliwe; 

3. W przypadku, gdy uzyskane w oparciu o metodę najmniejszych 
kwadratów oceny parametru a, procesu AR(1) charakteryzują się bez- 
względną wartością większą lub równą 1, to proces należy zaliczyć do 
klasy niestacjonarnych; 

4. W przypadku, gdy bezwzględna wartość oceny parametru aņ jest 
mniejsza od 0,925, to proces można zaliczyć do klasy stacjonarnych; 

5. Jeżeli bezwzględna wartość oceny parametru ay, należy do przedzia- 
łu (0,925; 1) nie ma możliwości określenia, do jakiej klasy należy proces. 
W takim przypadku należy proces ten rozważać jako proces niestacjo- 
narny. 
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PE3IOME 


B cTraTbe OGcyXKĄaIoTcA HpMMEHAEMBIE B HACTOAUIee BDEMA METOĄbI KCCJIEĄOBA- 
HMA CTAIIMOHAPHOCTM ONHOMEpDHBIX CJyUaAHBIX mpoyeccos, T.e. METOĄ Boxca- 
-Jlrkeaknica u MEeĄKAHHBIŃ TecT Cepuu. ONKPaACE Ha PEJYAbLTATBI CKMYJALWNKOHHBIX 
KCCIE4XOBAHHŃ, aBTOp OOGHAapyKMBAET  HECOBEPLIEHCTBO  HAJBAHHbIX MeTOJOB. 
B cTaTbe npen/laraeTcA HpoqeNypa, KOTOPAA NHpeĄCTaBIAETCA KOHKypeHTHOM TMO 
OTIŁOLIEHMIO K Ha3BaHHbIM BbiliO. IIpehjiaraeMaa npoqełypa npeącTaBiAeT coGoń 
PE3YJIbTAT CKMYJALNMOHHBIX MCCJIEXOBaHi, MpOBeAeHHbIX JAA umpokoro KJIACCA 
CTaAUMOHAPHbIX M HECTANUOHADHbIX ABTODETpECCMBHbIX IIDOLNECCOB. 
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SUMMARY 


The paper presents an evaluation of the methods which so far have been ap- 
plied to study stationarity of one-dimensional Stochastic processes. These methods 
are Box-Jenkins’ method and a median series test. The results of simulation stu- 
dies prove the deceptiveness of the methods under discussion. A procedure was put 
forward which seems to be competivite in relation to the ones mentioned above. 
This procedure is a result of simulation studies conducted for a wide class of sta- 
tionary and non-stationary auto-regressive processes. 


